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The PHOBOS experiment has measured the charged particle multiplicity at mid- rapidity in 
Au+Au collisions at -^/sjvjv = 200 GeV as a function of the collision centrality. Results on 
dN c h/dri\\ rl ^ <1 divided by the number of participating nucleon pairs {N part ) /2 are presented as 
a function of {N par t). As was found from similar data at -y/sjvjv = 130 GeV, the data can be equally 
well described by parton saturation models and two-component fits which include contributions that 
scale as N par t and the number of binary collisions N co u . We compare the data at the two energies 
by means of the ratio -R200/130 of the charged particle multiplicity for the two different energies as 
a function of (N part ). For events with {N part ) > 100, we find that this ratio is consistent with a 
constant value of 1.14 ± 0.01(stat.) ± 0.05(syst.). 
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Collisions of gold nuclei at the Relativistic Heavy-Ion 
Collider (RHIC) at ^/snn = 200 GeV offer a means to 
study strongly-interacting matter at high densities and 
temperatures. The goal is to create a large- volume, long- 
lived state within which quarks and gluons are no longer 
confined within hadrons, the quark-gluon plasma (QGP). 
The role of the collision geometry in determining the ini- 
tial parton configuration is important for understanding 
any collective effects which may be present in such col- 
lisions. We can study this by means of the mid-rapidity 
charged particle multiplicity as a function of the number 
of nucleons that participate in the collision, N part . Mea- 
surements of proton-nucleus reactions at lower energies 
suggested that the charged multiplicity from soft pro- 
duction mechanisms should simply scale with N part 
With increasing energy, however, one might expect some 
component of particle production to depend on the num- 
ber of binary collisions, due to the increasing role of hard 
processes (minijet and jet production). In nuclear colli- 
sions, (Ncou) oc (Npart) 4 / 3 , making these systems quite 
suitable for studying the interplay between the various 
effects. 



Nuclear collisions at RHIC also provide an opportunity 
to study Quantum Chromodynamics (QCD) in a novel 
regime where parton densities are high, yet the strong 
coupling constant is small due to the large momentum 
transfers involved. In such a regime, gluon densities can 
be large enough that the gluons recombine, causing a sat- 
uration of the gluon structure function at low Bjorken x, 
characterized by a momentum scale Q s . Since the parton 
densities in the initial state can be related to the density 
of produced hadrons in the final state, definite predic- 
tions are possible regarding the multiplicity of charged 
particles as a function of energy and centrality || . 

A recent extension of the calculations by Kharzeev and 
Levin Q has given predictions of the energy, rapidity, and 
centrality dependence of the charged particle multiplicity. 
These new calculations use the predicted QCD evolution 
of measured results, incorporating parameters extracted 
from inclusive deep-inelastic electron-proton scattering 
data H . Of primary importance in this treatment is the 
exponent A, which parameterizes the energy dependence 
of the saturation scale as Q 2 S oc (y/s) x . Kharzeev and 
Levin use this to predict that the energy dependence for 
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dN c hldr) will also scale as {\fs) x at high energies. Fur- 
thermore, they predict that the higher energy collisions 
allow events with larger impact parameter to be in the 
saturation regime. This affects the multiplicity from pe- 
ripheral events more than for central events, which al- 
ready have sufficient parton density at lower energies. 
Thus, they predict the multiplicity in peripheral events 
to rise slightly faster with energy than for central events. 

This effect should be contrasted with the much sim- 
pler two-component parameterization constructed to in- 
terpolate between proton-(anti)proton {j>p,pp) and cen- 
tral nucleus-nucleus collisions by incorporating contribu- 
tions which scale with the number of wounded nucleons 
as well as the number of binary collisions. The pseu- 
dorapidity density at mid-rapidity has been measured 
to rise by 14±5% between 130 and 200 GeV in central 
Au + Au collisions while interpolations based on UA5 
data on pp collisions at similar energies suggest only an 
8% increase for elementary collisions. Therefore, the ra- 
tio #200/130 = dN ch /dr]\2oo/dN ch /dr]\i 30 should decrease 
with increasing impact parameter. 

The predictions from the saturation and two- 
component scenarios have been found to be nearly indis- 
tinguishable as a function of centrality at ^snn = 130 
GeV and agree well with the published RHIC data 
from PHENIX @ and PHOBOS @Q. The authors 
of Ref. H find this to be a nontrivial consequence of the 
saturation formalism, although this ambiguity is unlikely 
to be reduced on the basis of the multiplicity data alone. 

We have performed a measurement of the centrality 
dependence of the mid-rapidity charged particle multi- 
plicity in Au+Au collisions at ^/s^n — 200 GeV using 
the PHOBOS detector and derived #200/130 as a func- 
tion of the number of participants. PHOBOS has previ- 
ously published results for the centrality dependence of 
dN ch /dr]/(~{N part )) at ^/sWn = 130 GeV Jh} as well as 
#200/130 for the 6% most central events ({N part } ~ 343) 
[f| . The present results are an extension to the previous 
measurements, using similar methods of analysis. 

The collision centrality is determined using the signals 
measured in two sets of 16 paddle counters (PP and PN) 
located at \z\ = 3.21 m with respect to the nominal inter- 
action point. For events at z = 0, these detectors mea- 
sure charged particles produced into 3 < \rj\ < 4.5. As 
discussed in Ref. [jl0| , we rely on the monotonicity of the 
paddle signal with the number of participants (verified by 
correlations with the PHOBOS zero-degree calorimeters) 
to extract (N part ) for a chosen fraction of the total cross 
section, based on the Glauber calculation used by the 
HIJING model f[^| . The dominant systematic error on 
this quantity reflects the uncertainty on the estimation of 
the total cross section observed by the PHOBOS trigger 
counters. By a study of events with a low number of hit 
paddles using a full simulation of HIJING events, we have 
estimated the efficiency for triggering on the Au+Au to- 
tal inelastic cross section to be 96 ± 3% for ^/JWn = 130 



GeV and 97±3% for = 200 GeV. The slight differ- 

ence between the two energies stems mainly from the in- 
crease in the width of the pseudorapidity distribution || . 
In order to reduce the error on the ratio of the multiplic- 
ity as a function of centrality, we have analyzed the simu- 
lations at both energies in an identical fashion. This has 
led to a slightly different efficiency at ^snn = 130 GeV 
than the one presented in Ref. [fL0|, which was 97 ± 3%. 

The charged particle multiplicity has been measured 
independently for each energy using the same tech- 
nique described in [flo| , based on counting 3-point tracks 
( "tracklets" ) in the spectrometer and vertex detectors. 
This approach permits some level of background rejec- 
tion relative to simply counting detector hits, as was 
done in jy]. We use a simulation based on HIJING 
events and a full GEANT simulation to study the effects 
of occupancy, combinatorial background and experimen- 
tal backgrounds. Due to the better granularity in both 
pseudorapidity and azimuthal angle in the spectrometer, 
its associated systematic error is 4.5%, to be compared 
with 7.5% in the vertex detector. As the particle density 
does not increase dramatically between the two energies, 
the systematic error is the same for the two data sets. 

In forming the ratio #200/1301 we have analyzed the 
data in each subdetector with the same correction fac- 
tors for both energies. This reduces the importance of 
the precision of the Monte Carlo simulations used to de- 
rive the corrections at each energy, since they directly 
cancel in the final value of #200/130- This is justified by 
the fact that the particle density does not change sub- 
stantially between the two energies. 

We present two forms of #200/130- The first is done 
for a chosen fraction of the total cross section, which 
requires no other input other than the trigger efficiency. 
The other is the ratio for a constant value of N part , which 
requires a model of the nuclear geometry. In both cases, 
a precise measurement of the #200/130 is achieved by the 
average of separate measurements of the ratio in each 
sub-detector. 

In Fig.|l|(a) we show the ratio of dN c h/ dri^^i at 
\J s nn — 130 GeV and 200 GeV for the same percentile 
of the total cross section, #200/130- The values in each 
bin are listed in the last column of Table |. This is found 
to be approximately constant for the 50% most central 
collisions. The grey band indicates the magnitude of the 
systematic error, which is symmetric around the shown 
data points. It reflects the uncertainty in the estimate of 
the relative uncertainty of the trigger efficiency between 
the two center-of-mass energies. 

The ratio of the charged particle multiplicity between 
^/sJfN — 200 and 130 GeV as a function of the number of 
participants has been obtained simply by taking the ratio 
of dN ch / 'c?77 1 1 T7 1 < x / '(^(Npart)) for the two energies, which 
corrects for the increase of the inelastic cross section be- 
tween y'SjvAr = 130 and 200 GeV (as parameterized in 
HIJING). The increasing cross section causes the ratio of 
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(Npart) between the two energies in each centrality bin, 
Npart/Nparti to increase by approximately 3% between 
the most central events and the bin corresponding to 45- 
50% of the total cross section ((N part ) = 65), as shown 
in Fig. |(b). 

It should be emphasized that in both ratio mea- 
surements, averaging over independently-measured ratios 
does not give precisely the same answer as that of taking 
the ratios of the averaged quantity at each beam energy. 
We have chosen the former method since it allows the 
correction factors for each measurement to cancel in the 
ratio. 

In Fig. |^(a) we present the separate results for 
dN ch /dri/{l (N part }) vs. Npart for = 200 GeV 

(shown as closed triangles) as well as 130 GeV (open tri- 
angles). For both energies, we observe a continuous rise 
of particle density with increasing centrality. The total 
95% C.L. error for the 200 GeV data is shown as a shaded 
band and includes both statistical and systematic errors. 
The errors at 130 GeV are not shown, but are of similar 
magnitude. InTable|[ we list (N part ), dN c h/ drj\\n\<i and 
dN c h/ dri\\n\<i/ {\(Npart)) for each energy and centrality 
bin. 

To put these results in context, we also show two cal- 
culations. The first calculation shown is the prediction of 
the parton saturation model for both energies, indicated 
by solid lines in Fig. |2|(a). The only parameter needed to 
predict the energy density is the exponent A, which is set 
to 0.25 as is done in Ref. 0. The calculation is truncated 
at (Npart) ~ 65 since this model is not appropriate below 
this value S. 

The second calculation involves fits to the data using 
the two-component parameterization proposed in Ref. 
@, dN ch /d V = n pp ((l - x)(N part )/2 + x(N coU )). The 
parameters have a simple interpretation as the fraction 
of production from hard processes (x) and the number 
of particles associated with a single pp interaction (n pp ) 
when Npart = 2 and N co u = 1. We have performed 
fits using a parameterization of (N co u) based on Ref. j|, 
(N co ii) = 0.352 x (Npart) 1 ' 37 . Using the measured val- 
ues for n pp we find that values of x — 0.09 ± 0.02 for 
= 130 GeV and x = 0.11 ± 0.02 for ^/IWn = 200 
GeV account well for the measured centrality dependence 
of dN c h/drj/(^ (Npart))- These are shown as dashed lines 
in Fig. |2](a). The error estimate is based on the al- 
lowed range of x which keeps the results of the calcu- 
lation within our stated error bands. The value of x 
obtained for ^/Jnn = 130 GeV agrees with the value 
extracted in Ref. ||, which only used the pp value and 
the central PHOBOS result from Ref. |§. The value 
of x at y/sMN — 200 GeV is consistent with the pre- 
diction in Ref. [Q that x should vary as the square 
of the gluon structure function, (y / SAT/v) 2A , which gives 
x(^Wn = 200) = 1.2 x x(^/sW^ = 130) for A = .25. 
These results are also in good agreement with the two- 



component calculations by Li and Wang |14[ ], which use 
recent parameterizations of the gluon structure functions 
as well as nuclear shadowing. 

The centrality dependence of -R200/130 is shown in Fig. 
||(b) and shown in Table Q. The error bars indicate the 
contribution from counting statistics, which becomes im- 
portant when the systematic effects cancel. The sys- 
tematic error on the centrality dependence of the ratio, 
whose magnitude is shown by the grey band, is symmet- 
ric around the shown data points and incorporates two 
major contributions. The first is the uncertainty in the 
relative trigger efficiency between the two energies, which 
is shown in Fig. |l|(a). The second error source is the 
change in N part for a given centrality bin as a function of 
energy, shown in Fig. 0(b). We have added half of this 
difference in quadrature with the first contribution to get 
the final systematic error. In addition, there is an overall 
scale uncertainty on -R200/130 and ^200/130 °^ which 
was discussed in Ref. || and stems from the Monte Carlo 
studies from which we derived our acceptance, feed-down, 
and background corrections. If we fit the ratio vs. N part 
to a constant -R200/130 above (N par t) — 100, we find that 
it is equal to 1.14 ± O.Ol(stat). 

For Np ar t — 2 we show the ratio of data from pp for 
yj snn = 200 GeV and an interpolation for y / sjv r /v = 130 
GeV 0. We also show the ratios of the saturation 
model predictions (solid line) and the two-component 
fits (dashed line) from Fig. ||(a). While the two cal- 
culations evolve in opposite directions as the impact pa- 
rameter increases, they remain sufficiently close down to 
(Npart) ~ 65 such that our present data cannot resolve 
them definitively. 

In conclusion, the PHOBOS collaboration has mea- 
sured the pseudorapidity density of charged particles pro- 
duced at mid-rapidity in Au+Au collisions at ^/snn = 
200 GeV. These data have been compared to similar 
data taken at ^Jsnn = 130 GeV by taking the ratio 
of multiplicities -R200/130 at a fixed value of (N par t)- For 
(Npart) > 100, we find that this ratio is approximately 
constant at 1.14 ± 0.01(stai.) ± 0.05(syst.), which is con- 
sistent both with the parton saturation model predictions 
and two-component fits to the data. 
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FIG. 1. a.) i?200/i30 vs - the average percentile of cross 
section. The grey band indicates the systematic error induced 
by the uncertainty in the relative trigger efficiency between 
the two beam energies, b.) The ratio N p ^ t /N p ^ t for each 
centrality bin derived using HIJING at -ysjvjv = 130 and 200 
GeV. 
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FIG. 2. a.) The measured scaled pseudorapidity density 
dN ch /dr]\ M<1 /(±(N part )) as a function of N par t for Au+Au 
collisions at ^/snn = 130 GeV (open triangles) and 200 GeV 
(closed triangles). The error band around the 200 GeV data 
combines the error on dNch/drjl^^i and {N par t)- The open 
and solid circles are pp results derived from the data in Ref. 
[7]. b.) The ratio of charged multiplicity for y^ijvjv = 130 
and 200 GeV, -R200/130! (f° r constant {N part }) vs. N part . The 
ratio from pp is shown by the open square. The grey band in- 
dicates the systematic error estimate. In both panels, results 
from a saturation model prediction and a two-component fit 
are shown as solid and dashed lines, respectively. 
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TABLE I. For each centrality bin, based on the percentile of the total cross section, we show the pseudorapidity density, the 
number of participants, and the scaled pseudorapidity density for both y / 5jvjv = 130 and 200 GeV. The errors shown on these 
values are systematic. We also show the ratio -R200/130 for each bin both corrected for the different (N par t) in each bin, and 
uncorrected ("raw"). The errors on the ratios are statistical only. 



Bin(%) 


dN ch /dri 


200 GeV 
(N par t} dN ch /dr)/{\ (N part )) 


dN ch 1 dr\ 


130 GeV 

(N par t} dN ch /dri/(\ (Npart)) 




Ratios 

-R200/130 


T>raw 
-"•200/130 





- 3 


700 ± 27 


358 ± 12 


3.91 


± 


0.20 


613 ± 24 


355 ± 12 


3.45 


± 


0.17 


1 


13 


± 


0.01 


1 


16 


± 


0.01 


3 


- 6 


629 ± 24 


331 ± 10 


3.81 


± 


0.18 


545 ± 21 


330 ± 10 


3.31 


± 


0.16 


1 


15 


± 


0.02 


1 


14 


± 


0.02 


6 


- 10 


548 ± 21 


298 ± 9 


3.68 


± 


0.18 


472 ± 18 


295 ± 9 


3.20 


± 


0.16 


1 


15 


± 


0.02 


1 


16 


± 


0.02 


10 


- 15 


455 ± 18 


256 ± 8 


3.56 


± 


0.18 


393 ± 15 


254 ± 8 


3.09 


± 


0.16 


1 


16 


± 


0.01 


1 


16 


± 


0.01 


15 


- 20 


376 ± 15 


217 ± 8 


3.48 


± 


0.18 


327 ± 13 


214 ± 8 


3.06 


± 


0.16 


1 


14 


± 


0.02 


1 


17 


± 


0.02 


20 


- 25 


312 ± 12 


183 ± 7 


3.41 


± 


0.18 


274 ± 11 


179 ± 7 


3.06 


± 


0.17 


1 


12 


± 


0.02 


1 


16 


± 


0.02 


25 


- 30 


252 ± 10 


152 ± 6 


3.32 


± 


0.19 


220 ± 8 


148 ± 6 


2.96 


± 


0.17 


1 


13 


± 


0.02 


1 


17 


± 


0.02 


30 


- 35 


202 ± 8 


124 ± 6 


3.25 


± 


0.19 


180 ± 7 


122 ± 6 


2.94 


± 


0.18 


1 


12 


± 


0.02 


1 


15 


± 


0.02 


35 


- 40 


164 ± 6 


103 ± 5 


3.19 


± 


0.21 


140 ± 5 


100 ± 5 


2.80 


± 


0.18 


1 


14 


± 


0.02 


1 


19 


± 


0.02 


40 


- 45 


130 ± 5 


83 ± 5 


3.14 


± 


0.22 


110 ± 4 


80 ± 5 


2.75 


± 


0.20 


1 


13 


± 


0.03 


1 


17 


± 


0.03 


45 


- 50 


95 ± 4 


65 ± 4 


2.90 




0.22 


83 ± 3 


63 ± 4 


2.64 




0.21 


1 


10 




0.03 


1 


14 




0.03 



•5 



